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Abstract 

General relativistic anisotropic fluid models whose fluid flow lines form a shear-free, irrotational, geodesic 
timelike congruence are examined. These models are of Petrov type D, and are assumed to have zero heat 
flux and an anisotropic stress tensor that possesses two distinct non-zero eigenvalues. Some general results 
concerning the form of the metric and the stress-tensor for these models are established. Furthermore, if the 
energy density and the isotropic pressure, as measured by a comoving observer, satisfy an equation of state 
of the form p = p{p), with ^ 7^ "jj then these spacetimes admit a foliation by spacelike hypersurfaces of 
constant Ricci scalar. In adclition, models for which both the energy density and the anisotropic pressures 
only depend on time are investigated; both spatially homogeneous and spatially inhomogencous models are 
found. A classiflcation of these models is undertaken. Also, a particular class of anisotropic fluid models 
which are simple generalizations of the homogeneous isotropic cosmological models is studied. 
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1. Introduction 



In a recent paper a comprehensive analysis of general relativistic spacetimes which admit a shear- 
free, irrotational and geodesic (SIG) timelike congruence was undertaken. This paper concentrates on a 
particular class of such spacetimes, namely, anisotropic fluid models whose fluid flow lines form a SIG 
timelike congruence (SIGA models). SIGA spacetimes are characterised by a single scale factor and a 
purely spatial three- metric If the additional assumption is made that the energy density and the 
isotropic pressure, as measured by a comoving observer, satisfy an equation of state of the form p = p[ii), 
with 3^ 7^ —5, then the Ricci scalar associated with the spatial three- metric is constant; we note that these 
spacetimes are not the same as the perfect fluid Friedman-Robertson- Walker (FRW) models, where the 
three-space is a space of constant curvature. However, the FRW models are a subclass of the SIGA models. 
The SIGA models studied in this paper are assumed to have zero heat flux and a non-zero anisotropic stress 
tensor, TTaf,, as measured by a comoving observer. The usual phenomenological assumption TTab — —Xaab 
does not hold for these models since the fluid congruence has zero shear. However, there are a variety of 
physical situations that are described by such SIGA models, for example, the interaction of a perfect fluid 
and a pure magnetic field 0, ||, ||, IJ]. Furthermore, since the anisotropic stress tensor is trace- free and 
three-dimensional in nature, it must either have three or two distinct non-zero eigenvalues. When there are 
only two distinct eigenvalues then the matter content of a SIGA spacetime is that of an anisotropic fluid; 
this paper is devoted to the study of such models. 

The paper is organised as follows: In section 2, we derive some general results about the form of the 
anisotropic stress tensor. The anisotropic stress tensor is shown to only depend on a single function, P, of 
the spatial coordinates, a three-dimensional vector that has unit norm with respect to the three-dimensional 
spatial metric and the spatial three-metric. The spacetime is also determined to be of Petrov type D. In 
sections 3 and 4, we specialise to the case where the three-dimensional hypersurfaces have constant Ricci 
scalar, ^R. In section 3, we use a 2-|-l split of the three-geometry to investigate a particular class of 
anisotropic fluid solutions which are simple generalizations of the FRW models. In section 4, we consider 
the class of solutions that have constant P and ^R, and prove that both spatially homogeneous and spatially 
inhomogeneous solutions exist. We finish by deriving some specific spatially homogeneous solutions for the 
case where both P and ^R are constant. 



2. The Stress-Energy Tensor and the Metric 

We consider general relativistic spacetimes with the following anisotropic fluid stress-energy tensor^ : 

Tab = ^iUaUb + PWUaUb + p^{UaUb ~ UaUb + gab) , (2.1) 

where m° is a unit timelike vector and n" is a unit spacelike vector orthogonal to Uan°- — 0. The scalars 
p\\ and p_L are the pressures parallel and perpendicular to n", respectively, and /x is the energy-density as 
measured by an observer moving with four- velocity m". Decomposing (2.1) with respect to ^"[^ we 



find that the energy flux relative to is zero, and that the isotropic pressure and anisotropic stress tensor 
(as measured by it°) are given by 

P = ^(P|| + 2p±) , (2.2) 

TTab = (P|| -P±){"a"b - -^(.gab + "a Mb)} • (2.3) 



*We shall follow the notation and conventions in Ellis [Q| and in in particular, all kinematical quantities are defined 
therein. Also, Latin indices range from to 3 and Greek indices from 1 to 3, and subscripts indicate differentiation with respect 
to the relevant spacetime coordinates. 
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In particular, when the timelike congruence formed by Ua is geodesic, shear-free, and twist-frcc then 
the line element of the spacetime may be written as ||^ 

ds^ = -dt^ + H'^{t)Kp{x'^)dx°'dx'^ . (2.4) 



The anisotropic stress-tensor ( p.3[ ) is now independent of t and is defined entirely in terms of the Ricci 
tensor, ^Rap, associated with the three- metric ha/3'. 



TTq/J = ^Rafi — -^ha/s^R ■ (2.5) 



Furthermore, tt^^ must satisfy the conservation equation 

y^'TTaf} = -^Vi3^R , (2.6) 
where is the covariant derivative with respect to the three-metric h^p- 

In terms of the coordinates (|]J) the vector Ua — {0,na) where Ua satisfies h'^^naUp — H^(t). 
Thus, the only non-zero components of the anisotropic stress tensor are 

7ra/3 = {p\\ -P±){nani3 - ^H'^{t)hai3{x'^)} . (2.7) 
The energy density fi and the isotropic pressure p are both functions of t only |^ , and are related to H by 

1 2 , 'i? _ + 'R 



/i 2H 6H^ + 3fl 



where 9{t) = 3H/H is the expansion of the timelike congruence. 

Theorem 1: If the stress-energy tensor has the form where the integral curves of form a shear-free, 
irrotational, geodesic timelike congruence then there exists a coordinate system in which the metric has the 
form (2.4). In addition, the non-zero components of the anisotropic stress tensor (in these coordinates) are 
given by (^.7[) where 



(i) p|| -p^ ^ Pix")/H\t) , (2.10) 
(n) na = H{t)Naix'^) , (2.11) 

such that Na is a unit vector with respect to the spatial three- metric hap. 

Proof: Without loss of generality the three- metric can be diagonalised ||. We take the metric ha/s to 
have the form 

haf3 = diag{A\B\C^) , (2.12) 

where A, B, C are independent of t. At least one component Ha is non-zero, say ni ^0. If n2 = (and/or 
"-3 = 0) then the component 7122 implies (p|| — p±)H^ = P{x'^) where P is as yet some undetermined 
function. If 77,2 7^ and rts 7^ then T^apix'^) — (p|| — p±)nanp when a ^ f3. Therefore, there exist functions 
fa^{x^) such that Ua/np — fa'^. Thus, the off-diagonal components of may be written as 

= iP\\-p±)inc.fU'' (a^/3) . (2.13) 

(There is no summation over the repeated index a in the above equation.) Hence, there exists non-zero 
functions Gq,(x^) such that 

(Pll -px)(ria)' - G„(x^) . (2.14) 
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Two possibilities now exist. Either tth 7^ or tth — 0. In both cases, (2.7) together with ( 2.12| ) 



and (2.14) implies there exists a function P — P(x'^) such that (p|| —p\^H^ = P. 



Part [ii) of the theorem now follows directly from part (z) and expression (2.7) for tTq^. 
Corollary 1: The anisotropic stress tensor may be written as 

7r„0 = P{x^){N^Np - hi^p} , 

where Na and P satisfy h°'^NaNp = 1, and 

iv„(3i?) = Vp{PN^)N^ - iv„P + PN^VpN^ □ 



□ 



(2.15) 



(2.16) 



The anisotropic stress tensor TTag has only two distinct eigenvalues which implies that the Weyl 
tensor for metric (2.4) is Petrov type D |]l^ |ll|. Conversely, if the metric is of the form (2.4) and tTq^ has 



two distinct eigenvalues then tTq^ must necessarily be of the form (2.15), and ( p.4[ ) can be interpreted as 
an anisotropic fluid, with four- velocity m° = 6f, whose flow lines are shear-free, twist- free, and geodesic. 
The proof is essentially as follows: let e"^-| be the two linearly independent eigenvectors corresponding to 
the same eigenvalue, A. Then the eigenvalue of the third eigenvector, e"-^^, is — 2A since 77^/3 is traceless. 
Furthermore, e"-^^^ is orthogonal to If e"j^^ is chosen to have unit norm then tTq^ can be decomposed as 



^a(3 _ 



(2.17) 



where g^^ gsc — 5^ and qab = G{A) • ^(b)- The desired result follows from the fact that the three-metric, 
hap, has the decomposition 

haj3 = e(i)Q 6(1)^ -I- g"*^ e(A)a e(s)^ . (2.18) 



3. A Particular Class of Solutions 



We wish to solve the three-dimensional fleld equations ( |2.5| ) and ( 2.15 ). As a general prescription, we 
write the field equations in terms of a 2-1-1 split. We suppose that we can choose coordinates {x,x^) such 
that iV^ cx dpx. (This will always be possible if is hypersurface orthogonal.) We define base vectors 
^Ia) ~ 9ax" and a lapse function M and a shift vector M'^ via 

— = MN- + Af^ef^) . (3.1) 

The extrinsic curvature and the intrinsic metric for the surface x =const are Kab and qab, respectively. 
They are related by 

d^gAB = 2MKab + '^^(aMb) ■ (3.2) 

It is unlikely that a general solution can be found. Hence, we impose some further physical constraints 
so that the matter content of the spacetimes under consideration is physically relevant. Specifically, we 
consider the case were the energy density and the isotropic pressure satisfy an equation of state of the form 
p = p{n)- Thus, either (i) =const or (ii) |^ = — ^ (see reference ||l|). We shall only consider the former 
case here, '^R — Rq — const. (For a discussion of the latter case see Q.) In addition, we shall only look for 
solutions that have unit lapse and zero shift. 
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Using the field equations ( pTs] ) and (2.15), and decomposing the Ricci tensor in terms of a 2+1 spht 
L2, nSl, we obtain the following equations: 



B = ^Rab + 2Ka^Kbc - KKab 
Ran - VbKa^ - OaK 

NN = -KabK^"" - d^K 



QxKai 



i^R - P)gAB/^ 
, 

(2P + 3i?)/3 , 



where K is the trace of the extrinsic curvature. The Ricci scalar is given by 



^R = ^R - - KabK^^ ~ 2d^K 



(3.3) 
(3.4) 
(3.5) 



(3.6) 



(Note that the Ricci tensor for any two metric is ^Rab — ^R9ab/2.) 

Under the conditions we have assumed here, equation (|3.2| ) can be integrated to yield 

<S>{x,x^)hAB{x^) 



9AB 



(3.7) 



with 2A = Equation (3^) implies that A = A(x) which in turn yields <& = 'ip{x) (l){x°'). We can 

absorb into Hab- Thus, without loss of generality, the two- metric (^^) can be expressed as 



9AB = i'{x)(l){x'-^)S_ 



AB 



(3.8) 



Now the function 2A = ip' (We have introduced the notation ip' = dxip-) For the remainder of the paper, 
we shall take — (y, z). 



The equation for the Ricci scalar (B.6) reduces to 

Ro = ^R - 6A^ - 4A' 



(3.9) 



The Ricci scalar for the metric (3.8) is easily calculated since the metric is conformally flat |l^, and is given 

by 

^^1-^ (3.10) 



'R 



(V^ is simply the flat space Laplacian). But ( |3.9| ) implies that ^R ~ ^R{x). Therefore, must satisfy 



V^ln^i + kef) ^ Q 



(3.11) 



where k is some constant. Solving (3.11) is equivalent to finding the conformal factor for a two-space of 
constant curvature. A rescaling of the y and z coordinates enables us to write the general solution as [|l6| 



[1 + -y + ^')]- 



(3.12) 



where the constant k has been normalised to take values ±1 and 0. 

Thus, inserting the results ( 3.12 ) and ( |3.10 ) into (3_^), we obtain the following differential equation 



for 



(3.13) 



We introduce a new function S{x) by letting ijj — S"^ . The function S then satisfies the differential equation 



2SS" + {S 



l\2 



(3.14) 
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The first integral of the above equation is 



2a i?o 2 



(3.15) 



where a is an integration constant. Equation ( 3.15| ) is very reminiscent of the Friedmann equation. Indeed 
an equation similar to ( |3.15 ) has been studied in the context of homogeneous and isotropic cosmological 
models p^, pi. 



An expression for P is found by subtracting (3^) from (3.5): 



5' 



52 



Equations ( ^.14 ) and ( 3.15 ) then imply 



3a 
S3 



(3.16) 



(3.17) 



which is in complete agreement with the conservation equation (2.16). 
Thus, the metric 



ds^ = ^dt^ + H\t) \ dx^ + S\x) 



[1 + 1(2/2 



^2)]2 



(3.18) 



can be interpreted as an anisotropic fluid, with four- velocity u°' — df, when k = ±1 or k = 0, and S{x) 
satisfies (3.15). The energy density, fj, and the isotropic pressure, p, as measured by a comoving observer 
are given by (2.8) and (|2.9|), respectively. The pressures pj_ and py are 

a 



P 
P 



m 53 

2a 

m 5*3 



(3.19) 
(3.20) 



We note that if a equals zero then the model does not represent an anisotropic fluid model since 
P — Q. This implies that py = pj_, and thus, the stress-energy tensor ( p.l[ ) corresponds to a perfect fluid 
and the spacetime is consequently a perfect fluid FRW model. Therefore, we shall only consider the case of 
non-zero a. 



A close examination of (3.15) reveals that it is remarkably similar to the zero pressure Friedmann 
equation. However, there is one notable exception; the constant a must be strictly non-negative for the 
analogue to be complete since a basically corresponds to the energy density, p, in the Friedmann equation. 
In our analysis, a simply represents an arbitrary integration constant, and thus may take on negative values. 



Equation (3.15) is difficult to analyse in general. However, a qualitative analysis can easily be 



undertaken. First, we shall consider solutions of ( 3.15 ) that allow S{x) to tend to zero as unphysical. In 



other words, there exists a number 5'min such that < S'min < S{x) for all values of x. In addition, if 
i?o > then S{x) must also be bounded from above, that is, there exists a number S'max < 00 such that 
S(x) < 5„ 



-'max- 



In terms of the Harrison- Robertson classification scheme ^7\, there are four basic types of 
solutions for which S{x) is non-zero everywhere. These solutions are referred to as: (i) static Si (unstable) 
and 5*2 (stable), where S{x) — Sg =const and S' = S" = 0; (m) asymptotic A2, Sg < S < 00 where 
Ss is a static point, that is, S'{S — Ss) — S"{S = Sg) = 0; (iii) monotonic M2, S'min < 5* < cxd with 
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S{x — ±cxd) — oo; and (iv) oscillatory O2, where S'min ^ S < S'max- The allowable solutions are summarised 
in Table 1. 



The static models are by far the simplest models. Equations (B.14) and (3.15) imply that S"^ = 
2k/Ro (provided _Ro 7^ and k ^ 0) and a — kRo/3. Thus, sign(i?o) = k- Equation (3. If) then gives 
P = — i?o/2 =const. Thus, the metric for the static models {k = ±1) is 



H^{t) \ dx 



dz^ 



[1 + f(y2+z2)]2 



(3.21) 



(The coordinates have been rescaled to normalise the Ricci scalar, Rq, to 2k.) The above metric is a 
Kantowski-Sachs metric [|l^ with equal scale factors. The energy density and the isotropic pressure are 
given by (2.8) and (p.9[), respectively, viz. 



3i/2 



P = 



6HH 



+ K 



3H^ 



The anisotropic pressures associated with the metric ( [3.21 ) are 

2HH + 



P\\ 



2HH + 



The following anisotropic fluid with planar symmetry is another particular example: 

ds"^ = -dt^ + H^{t){dx^ + S^{x){dy^ + dz^)} 



(3.22) 
(3.23) 

(3.24) 
(3.25) 

(3.26) 



and corresponds to taking n — Q. The only non-singular solution occurs when i?o < and a < 0, and then 

cosh <! W — *> . (3.27) 



S{x) = 



12a 
Ro 



4. Solutions with constant P 

In this section, we shall look for solutions that have both '^R and P constant. First, we note that this 



assumption greatly simplifies the conservation equations (2.16). However, if P is not constant then the 
resulting models are inhomogeneous; consequently a more pertinent reason for choosing this ansatz is the 
fact that we would like to determine the spatially homogeneous SIGA models. 

We shall assume that the matter associated with a spatially homogeneous model inherits its sym- 
metries |po| ; in particular, any scalars constructed from the stress-energy tensor should be independent of 
the spatial coordinates. For an anisotropic fluid this would imply that P is constant since P = P(x"). 
The vector Na must be both geodesic, N^W pNa = 0, and expansion-free, W^Na = 0, with respect to the 
three-metric hap (if P = const and "^R — Rq). We already know that there exist spatially homogeneous 
SIGA models with constant P and ■^i? (see [^). Thus, we now ask the following two questions:(l) are the 
Mimoso and Crawford solutions the only spatially homogeneous SIGA models, and (2) are there any 
spatially inhomogeneous models with P and '^R constant? 



7 



Thus, suppose that the Ricci tensor has the form 



(4.1) 



where the eigenvalues, a and /3, are constant and Na is a unit vector. The contracted Bianchi identities 
imply that Na can only have twist and shear, and hence 



(4.2) 



where ujap — '^[pN^] is the twist tensor and ggg — (pN^-^ is the shear tensor. If the shear tensor is 
zero then a theorem due to Bona and Coll implies that the three-metric hap admits a G4 as the 

maximal isometry group. If the shear is non-zero then we can use its unit orthogonal eigenvectors to form 
an orthogonal triad {A^q, -E(i)ai -E'(2)a}- The shear and twist can then be decomposed as 



E(2)a E(2)p\ , 
E{2)a £'(l)/3] , 



(4.3) 
(4.4) 



where the vorticity and shear scalars are given by 2ci;^ — uoafi^"'^ and 2a^ = cTapcr"'^ 7^ 0, respectively. 
The integrability conditions for Na, 



y^^lsNa - VpV^Na = N^ Rsc^fSj , 



can be employed to show that 



(-E'(I)q)' — (-E'(2)a)' 







(4.5) 



(4.6) 



where (• • •)' = N^Vai' ' ') the covariant derivative in the direction of N". Decomposing the covariant 
derivatives of i?(A)a in terms of the triad, we find 



V/3£'(i)q, 

V/3£'(2)q 



-6'2 i?(2)Q £'(l)/3 + ^1 -E(2)a -E'(2)/3 " O" £'(l)/3 " 1^ -^a -E(2)/3 , 
-9lE(^l)aE(2)l3 + 02E{i)aE(i)i3 + UJ E^i)^ + aNaE{2)l3 , 



(4.7) 



where Oi and 62 are the expansions of and £'(2), respectively. Equation ( |4.5| ) also yields three further 
equations: 



CTi 
0-2 



a = - 0-2) 

UJ2 + 2 a (6*1) = 
UJ2 + 2 a (6*2) - 



(4.9) 
(4.10) 
(4.11) 



where a\ — ^ Vq,(T, 02 — E%. V^a , etc 



Taking {iVo, dx", i?(i)c( dx", i?(2)a ^^2;"} as an orthogonal basis of 1-forms for the metric and using 
Cartan's structural equations, we obtain the extra integrability conditions: 



%)' + {^2)0 + (0i)cj = , 

9i)' - {Q2)0J - {Qx)CT = , 



92),2 + (^l).! + {^2? + (^l)' + /3 







(4.12) 
(4.13) 
(4.14) 



We now notice that ( 4.12 ) implies (Q2)" = —%{0 2)- In trodu cing coordinates = [x, y^) such that x is the 
affine distance along the integral curves of TV", ( 4.12 ) and ( 1.13 ) can then be integrated. At this stage the 
analysis splits into two parts: (i) a 7^ and (ii) a = 0. 



(i): If a 7^ then (4.14) can be employed to show that the only valid solution is 9i ~ 62 = 0, whence f3 — 0. 
Equations (4.9)-(4.11) then imply that both uj and cr are constants. Hence, the metric must admit a G3 



8 



1% since all the eigenvalues and spin coefficients are constant (see [E4j for a complete classification of 



all postive definite three-metrics that satisfy (4.1) and admit either a or a G4 as the maximal isometry 



group). Thus, we have proved the following result: 

Theorem 2: If hap is a positive definite Ricmannian three-metric whose Ricci tensor has exactly two 
distinct constant eigenvalues then the metric is homogeneous if either (i) the shear of the principal Ricci 
direction corresponding to the simple eigenvalue is zero or (ii) the degenerate eigenvalue is zero. □ 

We note that the Killing vectors of the three-space are also Killing vectors of the four-dimensional 
spacetime and hence, the resulting SIGA models are spatially homogeneous. 



(ii) If a — then we can always choose uj = a, and consequently ( 4.13 ) implies 82 — —9i = Q{y )■ Equations 



(|f.9D-(|4.1l|) yield 



0,1 - 0,2 - 2 + /3 , (4.15) 
CT 1 - CT,2 = 2 0-0 . (4.16) 

We note that N"' is not hypersurface orthogonal. However, we can pick a new orthogonal triad {iV", ^j^-j , 
where 



-(1) ^ ^(2) 



Cf±) = ± Ef^, /V2 , (4.17) 



such that ^"^^ is hypersurface orthogonal (we will label these surfaces by z = const) and such that 
is geodesic. If we let y be the affine distance along the integral curves of ^" -j then (x, y) will be intrinsic 
coordinates for the surfaces z = const and a 2-1-1 split can be employed to show that the three-metric can 
then be written as 

^ds^ = dx^ + dy^ + 2Midxdz + 2M2dydz + [M^ + {Mif + (Afa)^] dz^ , (4.18) 

where 

Ml = -y + g{z) and M2 = x + h{z) , (4.19) 

and 



A{z) eW-^y + B(z) e-V-l^y if (3 < 

M ^ { ^ ' ^ ^ (4.20) 

A{z) cos(V/32/) + B{z) sin(V^y) if /3 > 

where A, B, g and h are arbitrary functions. 

Theorem 3: If haf3 is a positive definite inhomogeneous Riemannian three- metric whose Ricci tensor has 
exactly two distinct constant eigenvalues then (1) the simple eigenvalue of the Ricci tensor is zero and (2) 



there exist coordinates {x,y, z) such that the inhomogeneous metric is given by (4.18)-(4.20). □ 

From theorems 2 and 3, we have the following corollary: 

Corollary 2: If ha/s is a positive definite Riemannian three- metric whose Ricci tensor has exactly two 
distinct constant eigenvalues then ha/s is not homogeneous if and only if both the simple eigenvalue of the 
Ricci tensor is zero and the shear of the principal Ricci direction corresponding to the simple eigenvalue is 
nonzero. □ 

We note that these spatially inhomogeneous models may be of some interest in the study of 



anisotropic fluid cosmologies. We also note that the three- metric, ( 4.18 ), possesses at most one Killing 
vector field (none in general) and that there exist examples of these spatially inhomogeneous models satis- 
fying all of the relevant energy conditions. 
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4..1 Spatially Homogeneous Solutions 



Let us explicitly construct some of the spatially homogeneous three-metrics. Now, since all of the spatially 
homogeneous cosmological models have at least one translational Killing vector field, say C = we shall 
consider the following class of three-metrics as an illustrative example: 

^ds^ = A^{y,z)dx^ + B^{y,z){dy^ + dz^) . (4.21) 



The Ricci tensor components Rxy and R^z for the above metric are identically zero, which implies 
that TTxy = TTxz = 0. Therefore, we must have P Ni N2 — P Ni N3 ~ 0. There are three possible solutions: 

(i) Ni^O and N2 ^ N3 = with P ^ 0, (ii) TVi = and P 0, and (in) P = 0. We can neglect the third 
possibility since it would imply that the stress-energy tensor represents a perfect fluid, as was noted earlier. 

jt) N,^0; N2 = Ns = 

The expansion is WNa = A^'^Ni^x = 0, and the acceleration equations are WNi = 0, V"A^2 = 
-N^A^^Ay = 0, V"iV3 = -NfA^^Az = 0. Hence, we can take Ni = A = 1 since Na is a unit vector. 
Thus, the metric belongs to the class of metrics discussed in section 3. Now, the metric must be the static 
metric solution (3.21) since the Ricci tensor has constant eigenvalues (both P and R constant). 

(ii) Ni = 

Now, A^i — ~ implies that there exists a coordinate transformation (y, z) {v, w) such that 
w is the afhn e-dist ance along the integral curves of A^"; that is, — 5^. Thus, since is geodesic the 
three-metric (4.21) may be written as 

^ds"^ = A'^{v,w)dx'^ + a^{v,w)dv^ + 2P{v)dvdw + dw^ . (4.22) 

The coordinate transformations w' = w + J I3{v)dv , v' — J P{v)dv reduce the above metric to the form 

ds^ = A^{v,w)dx^ + B^{v,w)dv^ + du? (4.23) 

(the primes have been dropped) and maintains the form of iV", A^" = 5'^. 

The expansion of A^Q, is V"A^c( = A^^A^ + B^^B^ = 0, which yields the equation 9tu(AP) = 0. 
Also, the field equations imply that R^^ = —A^^A^^ + A^^B^^A^B^ = 0, which implies that dwiA^/B) = 
0. Hence, we find 

= F{v) + G{w) , (4.24) 

r2 - g^*^") C4 25) 

' F{v) + G{w) ' ^^-^^^ 

where F, G, and g are arbitrary functions (we can set <? = 1 without loss of generality). Next, we can reduce 
the field equation R^^ = —B~^Bww ~ A^^ A^^ = (2P + Po)/3 to 

Therefore, we need only examine the two cases: (a) G = and (b) P = 0. 
(a) ^ = 



Here, A — A(v) — B^^ and equation (4.26) gives 2P + Rq — 0. Thus, B can be set equal to 1. The 
field equations R^x = {Ro — P)^^/3 and Rm — {Rq — P)B^/3 both reduce to the single ordinary differential 
equation 

A™ = -IRoA . (4.27) 
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We can rescale the {XjVjW) coordinates so that i?o — 2k. Then, the solution to (4.27) is 



. asin(ti) + 6cos(ti) if k = 1 

A = { ^ ' ^ ' , (4.28) 

aexp(w) + 6exp(— w) if k — —1 

where a and b are arbitrary constants. (When k = the metric is the flat FRW model.) If k = — 1 then 
there are essentially three distinct choices for A 

r (0 

A = < (a) cosh(u) . (4.29) 
I (iii) sinh(w) 

However, in each of the above cases, the two-surface w = constant is a surface of constant negative curvature. 
Thus the three metrics are diffeomorphic. The metric admits a G4 with a simply-transitive G3 subgroup 
of Bianchi type III and a G3 subgroup of Bianchi type VIII acting on the two-dimensional space w = 
constant. If k = 1 then the metric is a Kantowski-Sachs metric 

ds^ = -dt^ + H'^{t){sm^{v)dx^ + dv^ + dv?) . (4.30) 

Again, the three- metric is homogeneous and admits a G4. 

(b) F = 

In this case, A — A{w) = B^^. The field equations reduce to: 



(4.31) 



We observe that P — Ro, and that a solution only exists for _Ro < 0. Without loss of generality, we can 
take i?o = —2, and then we have the solution A = exp(±K7). The solution is the Bianchi type VIq metric 

ds"^ = -df + H^{t){e^'"dx^ + e-^^'dv^ + dw^) . (4.32) 

Mimoso and Crawford |2l] discussed Bianchi types I, III and Kantowski-Sachs metrics in the 
context of anisotropic fluid models. However, these are not the only SIGA models that can be interpreted 
as SIGA models. There are also Bianchi type II, VI, VIII and IX examples of SIGA models p4|. 





A 


+ A-^Al 


= (Ro- 


P)/3 


D-2 p _ 


A 




= {Ro- 


F)/3 


Rww ■ 




-2A-^Al 


= (2P4 


-i?o)/3 



5. Conclusion 

We have demonstrated that there exist coordinates such that the metric of a SIGA spacetime has the form 



{2A), where ^R^v{h-ys) = [a — P)N^Ni, + (3 /i^j,. In particular, we examined models in which both of the 
eigenvalues oi^Rap are constant; In this case, if either fi or the shear of is zero then the three-metric, h^v, 
is homogeneous and hence, the associated SIGA models are spatially homogeneous. A necessary condition 
for the three-metric to be inhomogeneous is a = 0. However, this condition is not sufficient. The static 
models discussed in section 3 are examples of spatially homogeneous models with a = 0. 

In this paper and our recent paper ||l[], we have only been interested in the study of spacetimes 
which admit a SIG timelike congruences. However, it is also of interest to study spacetimes which admit a 
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SIG congruence which is either null or spacelike. We shall investigate the null case in a future paper. In the 
spacelike case, perhaps the most natural (and most simple) setting in which such a congruence may take 
on a physical role is that in which the spacelike congruence is related to the preferred spacelike direction 
associated with an anisotropic fluid source [see ( ^.1[ )]. Let us now briefly describe this case since the analysis 
involved is similar to that employed in this paper. 



Thus, we study anisotropic fluid spacetimes with an energy- momentum tensor given by (2.1) in 



which the preferred spacelike vector n° orthogonal to is shear-free, irrotational and geodesic. In analogy 



with {2A), there exist coordinates in which n° — (5° and the metric is given by 

ds"^ = dx"^ + L'^{x)l^i:{x'f')dx^ dx'^ , (5.1) 

where ^,(,"0 = 0,2,3 [a;'^ = {t,y,z)] and l^ci^''') are the components of a Lorentzian three-metric with 
signature -|-1 (and coordinates can always be chosen such that is diagonal). Exact solutions of the 
Einstein field equations can now be sought in a similar manner to those obtained earlier in this paper, 
particularly in the case in which the adapted coordinates ( |3.1[ ) are also comoving coordinates (that is, in 
the case in which the anisotropic fluid four- velocity u"^ can be aligned along the timelike coordinate axis). 
However, we note that these exact solutions will be somewhat different in form to those obtained earlier 
due to the Lorentzian signature of l^^. We also note that in the special perfect fluid subcase these solutions 
are related to the solutions of Stephani |2^ . 

In particular, the results of Bona and Coll |2^ on Lorentzian three-metrics can be utilized in the 
special case in which the three-dimensional Ricci tensor, '^R^q, associated with the metric I^q has exactly 
two distinct constant eigenvalues. The results of section 4 can be generalized to the Lorentzian case and 
hypersurface homogeneous models can be obtained. 
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Table 1: 



K 


a 


Rq 


Type 





< 


< 


M2 


1 


< 


< 


M2 






> 


02,S2 


-1 


< 


< 


M2 




> 


< 


Si,A2,M2 



Table caption: Harrison-Robertson classification |g7|, |l8| for anisotropic fluid solutions with S'(a;) > 5min > 0. 
(The definitions of the constants k, a, Rq and the spatial types can be found in the text.) 



13 



